Results of wavelet processing of the 2i^'-capture ^^Kr experiment statistics 
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Results of a search for ^*Kr double if-capture with the large low-background proportional counter 
(2005-2008 years) at the Baksan Neutrino Observatory are presented. An experimental method 
and characteristics of detectors are described. Basic features of the digitized pulses processing 
using wavelet transform are considered. With due account taken of the analysis of individual noise 
characteristic it has been shown that the appropriate choice of both wavelet characteristics and 
sequence of processing algorithms allows one to decrease the background in the energy region of 
useful events with a unique set of characteristics by ~ 2000 times. New limit on the half-life of ^*Kr 
with regard to 2i4:-capture has been found: T1/2 > 2.4 ■ 10^^ yrs (90% C.L.). 
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I. INTRODUCTION 

Most /3/3 decay investigations have concentrated on the 
P~ f3~ transition (two-neutrino emitting and neutrinoless 
modes). This transition for the two- neutrino mode has 
been discovered in direct and geochemical experiments 
for ten nuclei of the "SCa, '^^Ge^^Se, ^^Zr, ^o^Mo, "^cd, 
i28Te, i30Te, 1^0^^ ^nd ^^^\J The data obtained 

for 2i^-mode offer a chance to directly compare different 
models of the nuclear structure, which form the basis for 
calculations of nuclear matrix elements |M^''|, and to se- 
lect the optimal one. Though direct correlation between 
the values of nuclear matrix elements for the two- neutrino 
and neutrinoless modes of (3/3 decay is absent lacking, the 
methods for calculating jAf^"! and |M'''^| are very close, 
and a chance possibility to estimate their accuracy in 
calculating |M'^''| appears only when comparing experi- 
mental data and theoretical results calculations for the 
probability of 2i//3/? decay. 

It can be expected that acquisition of experimental 
data on the other types of (3/3 transitions [the decay 
with emission of two positrons (/?"*" Z?"*"), capture of bound 
atomic electrons with emission of a positron {(3^ EC) and 
double capture of two bound atomic electrons {ECEC)] 
will make it possible to considerably increase the quality 
of calculations for both 2i/ and Ov processes. Much efforts 
have been currently made in searching for these processes 
[3,[1 and [i, in spite of the fact that the /3+/3+ and 
(3'^EC modes are strongly suppressed relative to (3~ /3~ 
decay due to the Coulomb barrier for positrons, and a 
substantially lower kinetic energy attainable in such tran- 
sitions. Positrons are absent in the final state of the 
2vECEC transition, and the kinetic energy of the tran- 
sition may be rather high (up to 2.8 MeV), which dictates 
determines an increased probability of a decay. However, 
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this process is also difficult to detect, since it is only char- 
acteristic radiation that is detectable in it. 

The experiment to search for 2_R'-caputre in ^*Kr using 
a copper low-background proportional counter of large 
volume has been carried out in the Baksan Neutrino Ob- 
servatory INR RAS since June 2005 |6l. The data ob- 
tained during first 159 hours of measurements with en- 
riched isotope of krypton allowed us to estimate the sen- 
sitivity {S) of the new installation in the experiment to 
search for ^^Kr 2Jsr-capture half-life. In one year of mea- 
surements it has reached S =1- 10^^ year (90% C.L.). 

Theoretical calculations based on different models give 
the following half-life times for this process in ^*Kr: 
3.7 • lO^i years 0; 4.7 • lO^^ years [|; 7.9 • 10^3 years 
Q. The last two values were obtained from the estima- 
tion of the half-life time in ^*Kr with regard to the total 
number of 2i/£^Ci?C-transition where 78.6% of events is 
due to 2i^2iir-capture [To| . Comparison of experimental 
and theoretical values shows that the sensitivity of the 
measurements has exceeded the lower limit of theoreti- 
cal calculations. This fact allows us to test the model 
of Aunola and Suchonen [7] . The technique to search for 
2iir-capture used in our research is based on the following 
considerations discussed below. 



II. BASIC ASSUMPTIONS 

When two electrons are captured from the if-shell in 
^*Kr, a daughter atom of ^^Se** is formed with two va- 
cancies in the K-sheW. The technique to search for this 
reaction is based on the assumption that the energies of 
characteristic photons and the probability that they will 
be emitted when a double vacancy is filled are the same 
as the respective values of the case when two single va- 
cancies of the if -shell in two singly ionized Se* atoms 
are filled. In such a case, the total registered energy is 
Si^ab = 25.3 keV, where Kah is the binding energy of a 
JT-electron in a Se atom (12.65 keV). The fluorescence 
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yield upon filling a single vacancy of the K-sheW in Se is 
0.596. The energies and relative intensities of the charac- 
teristic lines in the iiT-series are Kai — 11.22 keV (100%), 
Ka2 = 11-18 keV (52%), Kfn = 12.49 keV (21%), and 
K^2 = 12.65 keV (1%) The probability of deex- 

citation of a doubly ionized K-sheW through the emis- 
sion of Auger electrons (ca, Sa) only, or through a single 
characteristic quantum and an Auger electron {K, Ca), or 
through two characteristic X-rays and low-energy Auger 
electrons {K,K,ea) is pi — 0.163, P2 — 0.482 and 
P3 = 0.355, respectively. 

A characteristic photon can pass a long distance in the 
gas from the point of its origin to the point of its ab- 
sorption. For example, 10% of characteristic photons of 
12.6 keV energy are absorbed on a length of 2.42 mm 
{Pgas ~ 4.36 Bar) [l^]. Auger electrons of the same en- 
ergy will be absorbed on a length of 0.44 mm, which 
produce almost pointwise charge clusters of primary ion- 
ization in the gas [l3|- In the case of emission of two 
characteristic photons absorbed in a gas the energy re- 
lease will be distributed between three pointwise regions. 
It is these events that have a number of unique features 
and are the subject of study in this paper. 



III. EXPERIMENTAL TECHNIQUE 

To register this process a large proportional counter 
(LPC) with a casing of Ml-grade copper has been used. 
LPC is a cylinder with inner and outer diameters of 140 
and 150 mm. A gold-plated tungsten wire of 10 /im in 
diameter is stretched along the LPC axis and serves as 
an anode. The potential of -1-2400 V is applied to the 
wire, and the casing (the cathode) is grounded. Both 
ends of the anode are electrically connected to the to the 
corresponding end cap flanges via high-voltage pressure- 
sealed ceramic insulators with a central electrode taken 
from spark plugs. 

To reduce the influence of the edges on the operating 
characteristics of the counter, the end segments of the 
wire are passed through the copper tubes (3 mm in di- 
ameter and of 38.5 mm length) electrically connected to 
the anode. Gas amplification is absent on these segments, 
and charge is collected in an ionization mode. Taking into 
account teflon insulator dimensions, the distance from 
the operating region to the flange is 70 mm. 

The length of the LPC operating volume is 595 mm 
(the distance between the butt ends of the tubes), and 
the LPC operating volume is 9.159 1. The total capac- 
itance of the counter and outlet insulator is 30.6 pF. 
The total resistance of the anode and two output elec- 
trodes is ~ 600 Ohm. Indium wire is used to seal all 
detachable joints, and teflon gaskets are used to seal all 
nipple joints. The inner insulators are made of teflon, and 
their thickness was minimized to improve the degassing 
conditions during the vacuum treatment of the counter 
and to stabilize its operating characteristics. 

The LPC was placed inside the shielding of 18 cm 



thick copper, 15 cm thick lead, and 8 cm thick borated 
polyethylene layers. The installation is located in one of 
the chambers of the underground laboratory of the Gal- 
lium Germanium Neutrino Telescope experiment at the 
Baksan Neutrino Observatory, INR RAS, at a depth of 
4700 m.w.e. where cosmic ray flux is lowered by ~ 10^ 
times down to the level of (3.03 ± 0.10) x 10"^ cm^^g-i 

The counter is filled with pure krypton gas up to 4.42 
Bar having no quenching or accelerating gaseous addi- 
tions and is purified through a Ni/Si02 absorber from 
electronegative admixtures. Two samples of krypton 
were used in this work: one of 48.6 1 volume enriched 
in ^^Kr up to 99.8% [H and natural krypton (100 1) left 
after ^®Kr isotope extraction. Both samples have been 
specially cleared from cosmogeneous radioactive isotope 
^^Kr (Ti/2 — 10.756 yrs), present in atmospheric kryp- 
ton. 

The detector signals are passed from the anode wire 
to the charge-sensitive amplifier (CSA) directly attached 
to the high-voltage insulator. The CSA parameters 
have been optimized to transmit a signal with minimum 
distortions, and information on the primary-ionization 
charge spatial distribution in its projection to the counter 
radius is fully represented in the pulse shape. After am- 
plification in an auxiliary amplifier the pulses are col- 
lected by the LA-n20-12PCI digital oscilloscope. The 
oscilloscope, integrated with personal computer, records 
the pulse waveform as the numerical vector digitized with 
a frequency of 6.25 MHz. The length of the scanning 
frame is 1024 points (163.8 /xs) with a resolution of 160 
ns, ~ 50 /iS is "prehistory" and ^ 114 /is is "history". 

Processing of digitized pulses in " offline" mode is car- 
ried out using the technique developed in this work. 
This technique rejects pulses of non-ionzied nature, de- 
termines the type of event and its relative coordinate 
along the anode wire, and finally selects useful events 
according to chosen characteristics. 

IV. IDENTIFICATION OF MULTIPOINT 
EVENTS 

The response of a linear time-invariant system to some 
influence can be expressed by the convolution integral: 

oo 

g(t)= J h{T)l{t-T)dT, (1) 
— oo 

where in our case i{t) is the electric current at the out- 
put of LPC, q{t) is the signal from CSA, h{t) - the pulse 
transitive characteristic of CSA. To describe the output 
signal recorded by the digital oscilloscope in the form 
of a discreet set of instantaneous values of amplitudes 
one can move from integration to summation of instan- 
taneous function values with a step At: 

q{kAt) ^ At^h{nAt)i{kAt ~nAt), (2) 



' 50 55 60 65 

time, |is 

FIG. 1. A depiction of idealized pulse shape single-point event 
from LPC. The monotonically dashed curve a is a charge pulse 
pulse from the CSA output, and the front-back asymmetry 
solid curve b is is an electric current pulse (constructed by 
differentiating the charge pulse). 
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FIG. 2. A depiction of idealized pulse shape analysis parame- 
ters from three-point event expected from 2_R'-capture in case 
of registering two characteristic photon and a cascade of low 
energy Auger electrons from LPC: (a)- the electric current 
pulse, (6)- the charge pulse. 



The pulse shape of registered charged particles (or pho- 
tons) taken from LPC has some definite features which 
allow us to easily distinguish them from non-ionized 
noise signals and microbreakdown events. The registered 
events, in their turn, can be sorted out as "point-like" 
and "lengthy" events by comparing their pulse rise time. 
Point-like pulses with a small rise time for the pulse front 
are produced by particles which interact with a krypton 
gas in the counter in some particular region whose char- 
acteristic length is small in comparison with the cathode 
radius. Lengthy events with long of pulse rise time can 
be produced when several point-like events are simulta- 
neously registered at different distance from the anode 
("multi-point" or "multi-cluster" events). Besides, "ex- 
tended" pulses arise in a case when the fast electron track 
length is commensurable with counter radius. 

Idealized pulse shape events taken at the output of 
the LPC's and charge pulses taken from CSA's output 
produced by single-point events are shown in Fig.[T] The 
derivative of a pulse taken from CSA replicates the shape 
of a electric current signal from the LPC's anode wire. 
In a simplified case one can calculate the value of current 
in LPC by: 

^n = A, . (3) 

where Ax = x„ — x„_i. 

A pulse taken from the anode wire is formed mainly 
due to the negative charge induced onto the anode by 
positive ions moving towards the cathode and produced 
near the wire during the gas amplification (ionic compo- 
nent - i.e.). The contribution of the charge induced at 
the anode by electrons from avalanches (electronic com- 
ponent - e.c.) is small enough because they have to pass 
a small potential difference on their way to the anode. 
The electric current pulse shape is affected by the distri- 
bution in density of primary ionization electrons crossing 



the border of gas amplification region. The distribution 
parameters depend on drift time of the originally point- 
like charge towards the anode. The point-like ionization 
is smeared out during the drift time due to diffusion of 
the electrons into the cloud whose charge density distri- 
bution in its projection onto the radius is close to the 
Gaussian one. 

As is seen in Fig. [U the electric current signal pro- 
duced as a result of gas amplification has an asymmetric 
form and cannot be described by a single Gaussian curve. 
There is a mathematical procedure which allows conver- 
sion from a registered charge pulse to a electric current 
pulse produced by primary ionization electrons through 
the conventional border of gas amplification. The de- 
tailed description of the algorithm for electronic compo- 
nent segregation is given in [l6l |. The obtained form of 
a signal can be described by a set of Gaussian curves, 
thus determining the charge that was deposited in sep- 
arate components of a multipoint event. The calculated 
area of an individual Gaussian should correspond to the 
charge (energy) of the corresponding point-like ioniza- 
tion. Fig. [2] illustrates a demo version of a shape of a 
electric current pulse i{t): (a) from LPC and the corre- 
sponding CSA pulse ( b) from three-point event expected 
from 2iir-capture in case of registering two characteris- 
tic X-rays and a cascade of low energy Auger electrons. 
In Fig. [21 in addition to the primary pulse, there is a 
secondary pulse, smaller in amplitude, due to the pho- 
toelectrons knocked out of the copper case by photons 
produced in the electronic avalanches in the gas amplifi- 
cation process. 

Photoeffect on the cathode is probable enough because 
of absence of quenching additions in the krypton gas. The 
delay between pulse and after-pulse is determined by the 
total drift time of the electrons to move from the cath- 
ode to the anode. It presets the duration of the time 
interval to allocate totally any single event regardless of 
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FIG. 3. (a) and (6) curves demonstrate the converted pulses of 
voltage (charge) normalized for Ml of an original CSA pulse 
and obtained by integrating electric current pulses due to pri- 
mary ionization electrons (c) and (d) from two real events 
which are candidates for 2K-capture of ^*Kr. 



its primary ionization distribution over the LPC volume. 
In case of pure krypton the calculated drift time for the 
ionization electrons to move from the cathode to the an- 
ode is 58 fxs. A Gaussian area or a sum of Gaussians in 
case of a multipoint event for an interval of 58 fis, start- 
ing from the beginning of a pulse, gives total number of 
primary ionization electrons. 

For the pulse shape analysis we have used the following 
signal parameters indicated in Fig. O 1) T{ = tog — igi is 
the pulse rise time; ioi and io9 are the moments where 
the pulse amplitude has reached 10% and 90% of its max- 
imum value Ml; 2) Tp is the time of the after pulse ap- 
pearance since the beginning of the primary pulse; 3) 
A = {M2-M1)/M1 is the ratio of the difference between 
the after-pulse's and primary pulse's maximum ampli- 
tudes to the primary pulse amplitude. 

The parameters borders corresponding regions where 
useful signals are expected, and are determined from the 
distribution plotted for the pulses with known charac- 
teristics obtained from the calibration sources [31 . Us- 
ing these borders we select the candidate-events for 2K- 
capture from the whole set of events registered in the 
course of basic measurements. 

Fig. [3] (a) and (&) show converted pulses of voltage 
(charge) normalized for Ml of an original CSA pulse 
and obtained by integrating electric current pulses due 
to primary ionization electrons [(c) and (d)] from two 
real events which are candidates for 2iir-capture of ^^Kr. 

As is seen from Fig. |3] (c) and (d) the electric current 
pulses are very noisy. Noise and possible electrically- 
induced signals can both mask the low-energy component 
and produce a false one. Smoothing procedures applied 
to the initial pulse, with a following differentiation of it, 
result in deterioration of electric current pulse resolution. 
As to traditional frequency filtering methods with vari- 
ous window functions like filters of Hemming, Winner, 



Sawitsky-Goley and others [16| applied to the electric 
current signal, it does not allow one to always separate 
with enough reliability closely located other components 
of a compound event masking each other. 



A. Usifeul signal extraction from noisy LPC 
electric current pulses with wavelet analysis 

At present the mathematical base for wavelet analy- 
sis is well developed and present a good alternative to 
the Fourier transform in studying time (spatial) series 
with distinct irregularity. Wavelet transform, for a one- 
dimensional signal in particular, consists in its decom- 
position, through zooming and transfers, on a basis con- 
structed from a wave-like oscillation function with defi- 
nite properties called wavelet. Each function of this ba- 
sis is characterized both by a definite spatial (time) fre- 
quency and its localization in the physical space (time). 
The localized basis functions of wavelet transform resem- 
ble the signals under investigation to a greater degree 
than the Fourier basis functions. 

Thus, the distinction between the Fourier transform, 
conventionally used for signal analysis, and the wavelet 
transform is that the latter provides two-dimensional de- 
coupling of a regular signal where frequency and coordi- 
nate are considered to be independent variables, which 
means that in the wavelet transform an increase in res- 
olution of one of these variables does not automatically 
result in the decrease of resolution of the other. This 
independency of both variables allows one to simultane- 
ously analyze the properties of a signal both in physical 
(time, coordinate) and frequency spaces [13- [111- The 
wavelet transform has also a distinction from the short- 
time Fourier transform (the Gabor transform) resulting 
in a two-dimensional spectrum. The distinction is that 
one and the same function constitutes the basis of the 
former while in the case of the latter the functions are 
different [l7j . and thus the wavelet transform provides 
the better time-frequency localization, flexibility of anal- 
ysis and possibility to chose the best fitting form of the 
wavelet function. 

For our applications, the extraction of a real signal 
from noisy data is better approach mathematical tool of 
discrete wavelet transform (DWT). DWT algorithm is 
related to based upon uses the discreet vector {f{n)\n G 
N}, for which a wavelet spectrum is caluculated which 
in its turn is also a discreet vector. For the orthonormal 
wavelet bases on a binary lattice, there have been de- 
veloped algorithms of multiresolution analysis (MRA) or 
multiscale approximation (MSA). The idea of the latter 
consists in representing a signal by a sequence of images 
with different degree of detailing pL^]. Each image con- 
tains independent nonoverlapping information about a 
signal in the form of wavelet coefficients which are eas- 
ily calculated by the iteration procedure known as Fast 
Wavelet Transform. Taken as a whole, they solve the 
problem of the total analysis of a signal and substantially 
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simplify diagnostics of the underlying processes. 

MRA application allows one to represent the signal 
under investigation, x{n) , in the form of decomposition: 

m— J— jo m—0 

where ipj^^min) = \p2J°^(V°n — m) and '4'j.m{n) — 
V2^ip{2^n — m) are the known orthonormal scaling func- 
tions (pj^^rn{n) and wavelet functions ipj-min), respec- 
tively ; a,jg^m = {x,ipj„,m) are empirical approximation 
coefficients and dj^m = {S:,ipj.m) are empirical detail co- 
efficients; j^m £ Z are the current values of scale and 
shift; Nj^ (Nj) is the number of approximation coeffi- 
cients (detail coefficients) considered at the correspond- 
ing levels of decomposition; jo is the initial scale value; 
J is the final scale value. 

The resolution of the wavelet reconstruction x'^ {n) of 
an initial signal x{n) is set by parameter J. Eventually 
\\ x'^ — X >■ with J — ^ 00. 

Scaling functions {^j.m{n)} and mother wavelet func- 
tions {ipj,min-)}j^mez with 2^ non-zero coefficients satisfy 
the so-called two- level relations [n^ : 

25-1 

(p{n) = \/2 ^ h{k)ip{2n - k), 

k=0 
25-1 

k=0 

where h{k) and g{k) are the coefficients of low and 
high frequency filters of wavelet transform where g{k) = 
{-!)'' h{2^-k-l). 

The discrete orthogonal wavelet-transform at MRA 
can be defined as a recursive algorithm: 

0'{j-l),m = X] hk-2mO-jk, 
keN 

d{j-l)^m = (-'l)'^.9-fc-2m+lQjfc, 
k^N 

ao = X. 

According to the Mallat pyramid algorithm [20| , the ini- 
tial signal is first passed through the decomposition fil- 
ters of low and high frequency. The next step is to obtain 
approximation coefficients dj^m at the output of low fre- 
quency filter and detail coefficients dj^m at the output 
of high frequency filter, using decimation \, 2. This al- 
gorithm goes on according to the scheme given below 
until it reaches the decomposition level J. In terms of of 
wavelet coefficients the wavelet decomposition of a signal 
is presented in the following way. 

W{x) ^ dQ ^ {di.di} ^ {(22, d2,di} ^ ... 

... {dN, dN,dN-i, di} 

The solution of the noise reduction problem is carried 
out in four steps: 1) the initial signal is decomposed; 2) 



the threshold value for noise is chosen for each level of 
decomposition; 3) the threshold filtration of detail coef- 
ficients is carried out; 4) the signal is reconstructed. 

This technique is a non-parametric estimation of a re- 
gression model of a signal using an orthogonal basis (2l| - 
|23 | and works best for signals whose decomposition has a 
small amount of detail coefficients considerably different 
from zero. 

The reconstruction of a function x{n) is carried out 
using the modified wavelet coefficients for the signal de- 
convolution: 

^jo,m'Pja,m 

X! '^3.rnll^j,rn(n), (4) 
?n— j—jo rn—0 

where dj^^m- and dj^m-coefficients passed wavelet thresh- 
olding. 

As a rule, the criteria of entropy minimum (iJ) is used 
to choose the optimum wavelet decomposition, being a 
logarithm of a signal's energy H = J2 ^og (i^) . 

Entropy of an initial signal reaches its maximum due to 
noisiness of a signal. With the increase of the wavelet de- 
composition the entropy decreases to its minimum, which 
corresponds to the optimal level of wavelet decomposition 
of the initial signal. 

B. Determination of the optimal noise reduction 
for a signal from LPC 

When studying the effects constituting a small frac- 
tion of the background one should, apart from high qual- 
ity of experimental data, impose high requirements for 
their processing. Setting of the input parameters for the 
arbitrary noisy signals can present some difficulties in 
our case due to the necessity of having a priori data or 
an analysis of statistical characteristics of high frequency 
components of the signals. Since in the typical double- 
beta decay measurements (thousands of hours) the spec- 
trum of noise can be exposed to substantial changes 
due to possible electronic noise pickups and microphonic 
noise, some adaptive method of noise reduction is needed 
in order to minimize the uncertainty of the initial signal 
shape evaluation which substantially influences the de- 
termination of the type of event. 

It should be remarked that the choice of both a wavelet 
to be used and a level of decomposition depends on the 
properties of the signal under investigation. More smooth 
the wavelets create more smooth approximation of a sig- 
nal, and vice versa, "short" the wavelets trace peaks of 
approximated function is better. The level of decompo- 
sition affects the scale of sifted out details. With the in- 
crease of the level of decomposition the model subtracts 
the increasing level of the noise, with the result that pos- 
sible soothing occurs not only of the noise but also of 
some local characteristics of a signal. A series of approx- 
imation and detail wavelet coefficients is needed for the 
wavelet transform. 
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The noisy component of a signal in most cases is re- 
flected in the detail coefficients dj^m, and it is these co- 
efficients that are subject to processing in the noise re- 
duction technique. The noisy component, as a rule, has 
an absolute amplitude value less than that of the basic 
signal. So to reduce the noise one should zero those co- 
efficients that are less than some threshold value. The 
choice of the threshold level of noise affects the quality of 
the signal noise reduction, which could be evaluated as a 
signal to noise ratio (P): 



n=l 



n— 1 ^ 



(5) 



where s(n) and z(n) are the discrete vectors of a signal 
and noise; while iV^ and iV^ are their lengths, respec- 
tively. 

Setting a small threshold keeps a residual noise in the 
detail coefficients and results only in insignificant increase 
in the signal to noise ratio. With large enough threshold 
one can lose coefficients carrying essential information. 
Search for an optimum threshold amounts to finding the 
maximum value of the signal to noise ratio with the least 
quadratic deviation of the estimated signal. 

Let us express the empirical detail coefficients 



{(ij,m|m 
relation 



where 



1, A'j } corresponding to level j by the linear 



9], 



(6) 



J, ml 



l,Nj} are the true detail coefficients 
of a signal without noise {Nj in this case is the num- 
ber of detail coefficients, considered at a decomposition 
level j), and {'dj^m\^j,m G Norm(0, CT|),m — l,Nj} are 
the reading of the additive Gaussian noise with zero 
mean and variance, cr|. Then the solution of the noise 
reduction problem could be reduced to the search of 
estimates {gj,m\'m' — of true detail coefficients: 

gj,m = T(dj^m)- Such an estimate, carried out on the ba- 
lm 



sis of empirical coefficients {dj 



l,Nj} and a given 



threshold 9j^i, is in fact a construction of a regression 
model of the true coefficients: 



(7) 



where (j)(dj^rn,(^j,i) is the remainder term of thresholding 
function written in a general form. Expressing by the 
least-squares procedure the deviation of coefficients in 
the regression model ([7]) as some risk function 



Ki 



Hi 



hi) - X! ^3],k - 9],kf, 



(8) 



fc=i 



where i d N, we learn that the optimum value of the 
threshold 9j , producing in accordance with the Stein cri- 
terion [2^ the best noise reduction of a signal, corre- 



sponds to the case where function 
tremum: 

9j = arg _ min Rj{9j,i) 



has a global ex- 



(9) 



To meet this criterion used in the noise reduction we 
applied the Birge-Massart technique [IJ] to determine the 
threshold 9j to treat the detail wavelet coefficients. Coef- 
ficients dj^m less than the chosen threshold, were zeroed 
while the others were diminished by 9j. The optimum 
value of criterion 9j has been chosen in accordance with 
minimum " entropy - logarithm of signal energy" princi- 
ple. 

To enhance the reliability of both empirical signal 
shape determination and its parameters' extraction we 
can add to each initial frame of a noisy pulse a known 
model signal with preset parameters and then process it 
with wavelet analysis. As a rule, during first ~ 40 fis 
("prehistory") in the frame of digitized pulse from the 
CSA [see Fig. [3] (c) and (d)], there is a noise path only. 
In this interval of time we can select a vector of values 
containing no signal. Thus, we can produce a sampling 
control {z{n)\n — 1, Nx/4:;n = 1, iV2:/4} composed 
of two halves, each one being equal to "prehistory" {N^ 
is the total length of the frame). We add a model sig- 
nal So{n), composed of three Guassian curves summed 
together. We derive a model signal with a real noise hav- 
ing put in linear dependence of a width and center of the 
Gaussian curves on the rise time of CSA output pulse. 



s(n) = So{n) + z{n), n = 1, 



(10) 



To find the optimum technique for noise reduction of 
a similar signal and to define a quantitative measure of 
concordance between the estimated s{n) and initial So(n) 
signals, we introduce the following criteria: 



A 



N 

X; [s{n) ~ So{n) 

71=1 

N 

E So(") 

n=l 



N 

J2 s{n)so{n) 

B=^ ; (11) 

E sl{n) 

n=l 



C ^ ^A^ + {B -1)^. 



(12) 



The minimum value of functional A corresponds to the 
identity of the estimated and modeling signals, and its 
maximum corresponds to their total anticorrelation. In 
contrast, the minimum value of functional B indicates 
the anticorrelation of the mentioned signals and its max- 
imum standing for their total identity. So, the lesser A 
and greater B, the more accurate is the resurrection of a 
signal under investigation. In other words, the algorithm 
of the signal So(n) reconstruction amounts to minimiza- 
tion of functional A(s) with simultaneous maximization 
of functional B{s). The normalization of the introduced 
functionals allows comparing of various signals' decompo- 
sition for different wavelet function basis irrespective of 



the parameters of the signals and basic functions. Func- 
tional C(s) is a nonlinear combination of functional A(s) 
and functional B{s), and serves as a generalized quan- 
titative characteristic of a quality of the reconstructed 
signal and is to be minimized. 



C. Detection of a signal against additive noises 

The current task can be solved with the MATLAB us- 
ing Wavelet Toolbox [2^ for signal and image treatment. 
To choose the optimal wavelet allowing one to obtain the 
best of wavelet threshold of a signal, the criterion of max- 
imum ratio of initial entropy (Hq) and that purified from 
the noise signal (H) was used: 



V = Ho/H. 



(13) 



Confusing, you have to explain "tree" each type of 
wavelets using the criterion "entropy - logarithm of en- 
ergy" . The following orthogonal wavelets with com- 
pact carrier were used - Dobeshies (c?6N) , Simlet (symN) 
and Coiflet (coi/N), where N is the index number, 
designating the number of non-zero coeffeicients in fil- 
ters. Wavelets are not symmetrical and are not periodic 
enough. Simlet wavelets are close to symmetric to certain 
degree. 

For each specified empirical signal we defined the opti- 
mal number of decomposition levels using numerical sim- 
ulation for the case when the signal is resurrected to a 
large degree and the noise interference is not so large 
yet as to significantly distort the shape of a signal. The 
procedure of the signal resurrection is shown in Fig. |4] 
for the model signal Sq which has additively added noise 
interference Zo taken from a real signal x(n). As one 
would expect, both the quality of the estimated signal 
and noise-to-signal ratio depend on the degree of overlap- 
ping of the wavelet-spectra of a signal to be reconstructed 
and noisy signal. 

After choosing the best wavelet for a given pulse using 
the entropy criterion, the functional and (IT^ for 
each level of wavelet decomposition have been calculated 
for different values of threshold criterion of noise purifi- 
cation. Values of H, A, B and C are presented in Fig. |3] 
in (e), (/), {g) and {h), respectively, for 8 levels of pulse 
wavelet decomposition of s{n) [Fig. HTa) ] using wavelet 
of Simlet sym8 and soft thresholding [30j . 

Analysis of Fig. 0] (b), (c) and (d) allows one to ob- 
serve the change in the estimated signal with increase 
in number of decomposition levels L taken into account. 
Fig. ^c) shows that for decomposition level L — 4 oi a 
noise signal s{n) the minimum of functional C(s) corre- 
sponds to the maximum accuracy in model signal So{n). 
With the decomposition level of L = 3 we obtain the sig- 
nificant increase in sinusoidal noise, in case of (L — 5), 
which corresponds to the minimum of entropy H , the es- 
timated signal takes smoother shape resulting in stronger 
distortions of the initial model signal parameters. 
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FIG. 4. Signal modeling graphs with the added real noise (a) 
and de-noise of signals for various levels of decomposition: 
(fe)-L = 3, (c)-L = 4 and {d)-L = 5. (e) and (/), (g), (h) 
dependency's of entropy H and coefficients A, B, C from 
wavelet decomposition level, respectively. 



Information on the type of event is contained in the 
shape of a electric current signal. Therefore the mini- 
mum of functional C(s) approaches in our case the quan- 
titative criterion describing a relation between an initial 
modelling signal without noise and estimated signal is 
better. 

Thus, by determining, with the help of a model signal 
with definite empirical noise, the depth of the wavelet de- 
composition and threshold values 9j for each decomposi- 
tion level we find the resurrected empirical signal proper 
to the event in LPC. By describing signals purified from 
noise with the help of a set of Gaussian curves using the 
minimization root-mean-square error technique one can 
separate multipoint events from single-point ones. 



V. RESULTS OF THE FIRST STAGE OF 
MEASUREMENTS USING LPC 

Fig. [5ja) shows the total spectra of Ml amplitudes, 
normalized for 1000 h, obtained for the background of 
LPC filled with krypton enriched in ^^Kr (total collection 
time is 8400 h; spectrum " r , dark line) , and with natu- 
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FIG. 5. Pulse height Ml spectrum's (normalized for 1000 
h) of the background of LPC filled with krypton enriched in 
TSj^j. p ^ ^^^^ Yine) and with natural krypton ("^', light 
line): (a) - all events, (6) - single-point events, (c) - two-point 
events and (d) - three-point events. 



ral krypton (total collection time is 5000 ; spectrum " ^' , 
light line). Comparison of these spectra shows that the 
background of the counter filled with krypton enriched 
in ^^Kr exceeds the background with natural krypton at 
energy > 18 keV. The shape of difference (subtracted) 
spectrum in the range of 20 — 100 keV is described reason- 
ably well by the model spectrum composed of /3-spectra 
due to ®^Kr and ^''C decay in the proportion 1:17. Back- 
ground counting rates in this range are 56.3 h~^ and 47.2 
h^^ when filled with ''^Kr and "'^'Kr, respectively. The 
difference in ^^Kr isotope activities in the samples is due 
to its different residual value achieved during the isotopic 
purification. ^"^C isotope comes into the krypton gas dur- 
ing the continuous usage of the counter, supposedly as a 
result of slow sublimation of organic molecules (ethyl al- 
cohol and acetone were used to clean the components 
when assembling the detector) from the surface of the 
body of the counter. Overactivity of ^^C in the enriched 
krypton is due to the fact that this sample was mea- 
sured first and the main portion of this organic vapour 
penetrated into it. Actually, all the events in the operat- 
ing energy range produced by /3-particles absorption can 
be considered as single-point ones, except for those event 



where /3-particles have lost partially their energy through 
generating bremsstrahlung photon. 

At energies < 18 keV, the counting rate of events in 
spectrum is much higher than in spectrum "i". 
This is due to the presence, in the original atmo- 
spheric krypton, of cosmogeneous radioactive isotope 
^^Kr (Ti/2 = 2.29 • 10^ yrs) [Hi with volume activity of 
~ 0.1 min~^Z~^Kr 27|-[2^; the significant part of which 



comes into the sample of natural krypton during the pro- 
cess of its production. ^-'^Kr decays by electronic capture 
producing ®^Br (if -capture 87.5% [29]). i^T-capture yields 
a 13.5 keV energy release. One can see the peak of this 
line in spectrum (light line in Fig. [S]). The resolu- 
tion of the peak (FWHM) is 15.3%. Counting rate of 
events in (13.5 ±3.0) keV energy range is 220 h^^, which 
corresponds to the volume activity of ^^Kr (0.10 ± 0.01) 
min^^^^^Kr. 

The filling of if-shell vacancy of a daughter atom of 
bromine in 61.4% goes with omittance of characteristic X- 
rays of 11.92 keV {K^i, 100%), 11.88 keV (i^„2, 51.9%), 
13.29 keV {Kpi, 13.6%), 13.47 keV {Kp2, 1.36%) energies 
[26| and accompanying them Auger electrons of 1.55 keV, 
1.60 keV, 0.27 keV and 0.01 keV, respectively (relative 
intensities of i^Q^-lines are given in parentheses). Thus 
it is clear that in the case of photon emittance -ftTai, Ka2 
energy of Auger electrons is large enough to produce a 
distinct two-point event. Thirty-eight point six percent 
of all the events are due to the filling of Br atom Jf-shell 
vacancy accompanied by a cascade of Auger electrons, 
and are considered as events producing one point energy 
release. Events with the emittance of characteristic X- 
rays K^i and Ka_2 are also considered as one point events 
since energy release of Auger electrons do not exceed the 
CSA noise. Using the above data and taking into ac- 
count the efficiency of characteristic X-rays absorption 
in a krypton gas of LPC [Sp — 0.869) one can calcu- 
late the component analysis of complete absorption peak 
with energy 13.5 keV with 49.4% single-point -I- 50.6% 
two-point event [ld|. 

At 46.5 keV, in spectra "i" and [Fig. [^a)], one 
can see the peak corresponding to the source of ^^°Pb 
(Ti/2 = 22.2 yrs, ;3-decay, = 46.5 keV, with a yield 
is 4.25% per decay ^!|). ^^"Pb istope is produced af- 
ter 222Rn (Ti/2 = 3.82 days, a-decay) in the ^ssu de- 
cay chain. It can be produced directly in copper mate- 
rial due to trace radioactive contaminant decay (volume 
source) or accumulate on the surface of the counter body 
in the form of daughter isotopes of ^^^Rn decay, ^^^Rn be- 
ing present in the atmospheric air (surface source), dur- 
ing the preparation of mounting the detector. Gamma- 
radiation comes into a fiducial volume of the counter 
from both sources. There are mainly single-point events 
(due to photoeffect on krypton with de-excitation due 
to Auger-electrons) and two-point events (due to photo- 
effect with emission of krypton characteristic radiation) 
in the peak line of 46.5 keV, and only a small portion 
of events will be regarded as three-point ones where the 
primary photon or characteristic photon absorption oc- 
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curred through the scattering on the outer electrons with 
the following absorption of the secondary quantum. 

Conversion electrons (c.e.) can originate from the sur- 
face source, ^lOpb. 30.1 keV (52% per decay) and 43.3 
keV (13.6% per decay) lines are the most intensive. There 
are peaks in spectra 1 and ^ at '-^ 28 keV and ~ 41 keV 
that could be associated with these lines. All the events 
of these peaks should be single-point ones. The observed 
shift of maximum energy with respect to the expected 
line could be explained by the deposit of a portion of 
near- wall primary ionization electrons on the cathode due 
to their diffusion. With the emittance of conversion elec- 
trons the residual excitation of the ^^"^Bi daughter shell 
deenergizes by radiation of characteristic photons of L- 
series {El^ = 10.8 keV, 9.3% per decay; El^ = 13.0 
keV, 11.2% per decay and etc.) or/and by Auger elec- 
trons. In different combinations this radiation can enter 
the fiducial volume simultaneously with ^S-particles and 
c.e. In the simplest case, when c.e. with 30.2 keV and 
characteristic i^-photon are registered, a two point event 
with energy deposit of '-^ 42.5 keV is produced in the gas. 
The fraction of such events is small enough in the total 
spectrum. There is a step at ^ 34 keV [Fig. [5] (a) and 
(&)] corresponding to the escape peak of 46.5 keV line for 
krypton [E^ - Ek^^, = 46.5 - 12.6 = 33.9 keV). 

In general, all types of distortion of spectrum i, at the 
energies lower than ~ 50 keV, from the smooth drop- 
down base which is characteristic for energies higher than 
50 keV, could be "assigned" to ^^°Pb source. 

The scenario described above is confirmed by the dis- 
tribution of events related to different peaks in the ampli- 
tude spectra, composed separately for events with differ- 
ent number of point ionization clusters in the composite 
pulse. In Fig. [S] one can see the corresponding spectra of 
amplitude M\ plotted after selecting the type of events 
by analyzing the form of electric current pulses purified 
from noise: (6) - single-point, (c) - two-point, [g) - three- 
point. 

Comparison of the obtained relation of peak areas (13.5 
keV and 46.5 keV) for different components with esti- 
mated values allow one to determine the efficiency, £fc, 
of the procedure of the computer event selection with a 
given number of ionization regions. The efficiencies, e^, 
of two-point event selection were found to be 0.567 and 
0.733 respectively, which are in good agreement with cal- 
culated model values - 0.604 and 0.745 fl^. 

As seen from Fig. [SJ^), the main difference for back- 
ground spectra for LPC filled with different gases is 
in single-point (event) spectra. LPC background spec- 
tra with various gases are in good agreement with each 
other for two- and three-point events in the energy re- 
gion of 20 — 80 keV. For further analysis the three-point 
event spectra. Fig. ^d), have been used. Each event of 
these spectra is characterized by a set of energy deposits 
distributed over three point-like regions of the operat- 
ing volume of the counter. Energy deposits are propor- 
tional to the Gaussian areas (Ai, A2, A3), which describe 
three pointwise components of the total/composite cur- 
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FIG. 6. The three-point spectra selected under the conditions 
"Ci" and "Ca". 



rent pulse of primary ionization electrons on the border 
of the gas amplification region. Gaussian numbering cor- 
responds to the arrival time sequence of the components 
coming into the gas amplification region. To simplify the 
selection of events with a given set of features the ampli- 
tudes of components for each event are arranged in the 
increased order [(^1,^2, A3) — > (go < 91 < 92)]- 

In the sought-for events the minimal amplitude value 
(two groups of Auger-electrons from residual excitation 
with energies Ei ~ 1.5 keV and E2 ~ 2.9 keV) will corre- 
spond to a larger extent with respect to the resolution to 
the condition 0.9 keV < go < 4.5 keV: " CT . The mid- 
dle and highest amplitudes are produced by characteristic 

A'-photOnS {Kal + Kal,Kal + Ka2,Ka2 + Ka2, K^l + 

Ki3i,Kai+Ki32,Ka2+Ki3i,Ka2+Ki32) Sc**. In the total 
number of possible combinations, the fraction of the enu- 
merated combinations of the two photon (afe), registered 
in a composition of three-point event is ak = 0.985. It's 
convenient to introduce a parameter of average to max- 
imum amplitude ratio. With an allowance for resolution 
significant part of the event would have this parameter 
within 1.0 > qi/q2 > 0.7: " C2'\ The selection of events 
from spectra in Fig.[5I^d) corresponding to the conditions 
(" Cr'+" C2") allows an additional decrease in the back- 
ground of the expected peak energy region and gives a 
general representation of a spectra form for wide range 
of energies. 

The selected spectra are plotted in Fig. [6l They show 
that there are similar peaks for both spectra in the energy 
region of ^ 26 keV corresponding to the sought-for effect. 
To find their nature the study of the distribution of the 
number of background events of various types along the 
anode wire has been performed. This distribution was 
plotted using the dependence of the relative amplitude 
of the first after-pulse on the distance between the point 
of origin of the main pulse and the middle point of the 
length of the anode wire. This distance defines the solid 
angle viewing the inner surface of the copper cathode 
cylinder. The solid angle through which the operating 
surface of the cathode is viewed from the middle point of 
the anode wire and from the end points of the operating 
length of the anode wire are 3.9tt and ^ 27r, respec- 
tively [Fig.[71[a)]. The density of distribution of the solid 
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FIG. 7. (a) - the solid angle under which the interior surface 
of the counter from different points of an anode wire is visible. 
(b) - the density of distribution of the solid angle for the points 
uniformly distributed along the anode wire. 



angle for the points uniformly distributed along the an- 
ode wire is presented in Fig. [71[5). There is a dependence 
on the solid angle of the relative number of photoelec- 
trons, knocked out of the copper surface by the photons 
produced in the gas ionization. First after-pulse is pro- 
duced as a consequence of gas amplification of secondary 
photoelectrons generated on the cathode within the op- 
erating length of the anode wire. It is the A parameter, 
equal to the ratio of after-pulse and pulse amplitudes, 
(M2 — M1)/M1, accurate to the precision set by energy 
resolution of pulse and after-pulse, that determines the 
coordinate of the primary event with respect to the mid- 
dle point of the anode wire. To calibrate the counter with 
respect to A one needs a radioactive source uniformly dis- 
tributed over the LPC volume. ^^Kr is well suited for this 
task. 

A normalized distribution of the number of events of 
11.5 — 15.5 keV energy (^^Kr) from the spectrum of one 
point events of the LPC background with natural krypton 
versus parameter A is given in Fig.lHJa) as histogram with 
its maximum at A = 0.28. Distribution of the number of 
event with 20 — 80 keV energy range is given in the same 
figure in shaded region. 

Comparison of the distributions shows an excess of 
events over a bar graph in the region of A < 0.225 value, 
corresponding to the ends of the anode wire. The form 
of distribution of this excess is obtained by fitting the 
subtracted out histogram ["shaded region" minus "his- 
togram"], and is presented by dash-dot curve. Hence it 
follows that at the ends of the operating length of the an- 
ode wire there is an extra background source additional 
to the one that is uniformly distributed over the volume. 
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FIG. 8. The distributions of the events in LPC filled with nat- 
ural krypton with energies 11.5 — 15.5 keV (*^Kr, histogram) 
and with energies 20 — 80 keV (shaded region). The dash-dot 
and dotted curve is result of fitting operation of the shaded 
region, (a) - single-point events, (fe) - two-point events and 
(c) - three-point events. 



In Fig. |5] (6) and (c), there is similar distribution of the 
number of two-point and three-point events with energy 
20 — 80 keV (shaded region) and their uniform and near- 
boundary components (dotted and dash-dot curves). It 
is clear that the relative contribution of near-boundary 
into the corresponding spectrum increases with the in- 
crease of the number of point-like ionization clusters in a 
given event. 

The additional events at the ends of the operating 
length of the anode wire are produced by particles com- 
ing out of the 'dead' near-end volumes of a krypton gas. 
Apart from the proportionally amplified components a 
total pulse could also be composed of components of 
larger energy deposits that occur in the 'dead' volume 
and are collected in the ionization mode at the end bulbs 
of the anode. Specifically, the ~ 26 keV in spectra of 
Fig. m^a) could well be formed by alpha-particles emit- 
ted from the surface of the cathode in the near-end region 
where the anode wire is thickened with copper tubules. 
Alpha-particle can ionize if-shells of krypton atoms, and 
in this case the characteristic photons would have ener- 
gies of 12.6 keV. These photons can penetrate into the 
operating volume of LPC where gas amplification occurs. 
Registration of two photons gives two points. The third 
point is produced by ionization of alpha-particle collected 
in the ionization mode. Discarding pulses with A < 0.225 
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FIG. 9. The three-point spectra selected under the conditions 
{"Cr'+"C2") and A > 0.225. 



one can completely eliminate such events from 'three- 
point' spectra. The selection coefficient of useful event 
{k\ ) remains high enough: k\ = 0.84. 

In Fig. HI shaded-to-zero the bar graph (^^Kr) and the 
linear bar ("°*Kr) are three-point spectra selected under 
condition of A > 0.225 from spectra shown in Fig. [HI 
The background in energy region of (25.3 ± 3.0) keV in 
these spectra is ni = (38.1 ± 6.3) yr~^ and n2 — (20.7 ± 
6.0) yr""'^, respectively, which gives count rate of 2K- 
capture for ^^Kr to be Uexp = rii — n2 — (17.4 ± 8.7) 
yr~^. The obtained positive value does not exceed two 
standard deviations and in this case limit on the expected 
effect only was set for 90% C.L.: n^xp < 31.7 yr^^. 

The half-life limit has been calculated using formula 

Tl/2 > (^"2) • iV • P3 • Ep • £3 • Qffe • kx/Uexp, 

where N = 1.08 • lO^** is the number of ^^Kr atoms in 
the operating volume of the counter, pa = 0.355 is a 



portion of 2i4r-captures accompanied by the emission of 
two K X-rays; = 0.809 is the probability of two K- 
photon absorption in the operating volume; £3 = 0.422 is 
the selection efficiency for three-point events due to 2K- 
capture in ^^Kr; = 0.985 is the portion of events with 
two if-photon that could be registered as distinct three- 
point events; k\ — 0.840 is the useful event selection 
coefficient for a given threshold with respect to A. The 
result obtained is 

'Yi/2{'du + 2u,2K) > 2.4- lO^Vr (90% C.L.) 

VI. CONCLUSION 

We have demonstrated the possibilities for reducing 
the background of the proportional counter by ^ 2000 
times for the case of registering 2X-capture events in 
^*Kr by selecting pulses according to the number of point- 
like clusters and the event's coordinate along the length 
of the anode wire. In the spectrum of selected three- 
point events for LPC filled with enriched krypton, for the 
energy region of the sought-for effect, there is an excess 
of events that, within the achieved level of statistics, does 
not exceed two standard deviations. The current result 
has been used to set the limit on the half-life of ^®Kr 
2_Rr-capturc. The experiment is in progress. 
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